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Abstract 

The paper addresses rotor whirl stabilization of 
high performance turbomachinery which operates at 
supercritical speed. Basic whirl driving forces are 
reviewed. Stabilization and criteria are discussed. 
Damping seals are offered as a solution to whirl and 
high vibration problems. Concept, advantages, retro- 
fitting, and limits of damping seals are explained. 
Dynamic and leakage properties are shown to require a 
rough stator surface for stability and efficiency. 
Typical seal characteristics are given for the case of 
the high pressure oxidizer turbopump of the Space 
Shuttle. Ways of implementation and bearing load 
effects are discussed. Feasibility was first proven 
analytically in a Center Director Discretionary Fund 
project in fiscal year 1981. The Rocketdyne Division 
of Rockwell International has since installed such 
seals in the high pressure oxidizer turbopump for 
development to raise the speed limit and thus achieve 
higher thrust for the Space Shuttle Main Engine. 

Introduction 


High performance turbomachinery operates super- 
critically, i.e. with speeds above rotor resonances 
(critical speeds) . Thus more power is being packed 
into lightweight units unless speed must be limited to 
prevent a rotor from whirling and causing damage. The 
speed limit challenged many as well as the writer to 
find a solution, not without earnest prayer (which God 
answered,) , and not without the foundation layed by 
others such as Dr. Henry Black from Scotland, who is 
well remembered. *1-4 

* Deceased January 17, 1980. 
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Constant Whirl Forces 


Whirl is an orbital motion s superimposed upon the 
rotation of a rotor. Whirl forces are conceptually 
known s but accurate data are not readily available. 
Whirl forces (Fig. 1) are constant or increase with 
rotor motion. 

Synchronous constant forces are caused by mass 
unbalance, runout rubbing, and seal runout. Sub syn- 
chronous constant forces are generated by rotating 
Coulomb friction forces in rotor joints. 5 The con- 
stant forces produce a forward whirl limit cycle 
response at synchronous or sub synchronous frequencies. 
Well balanced rotors with minute runout yield low 
vibrations, however. Coulomb friction in rotor joints 
can cause damaging rotor vibrations. Coulomb friction 
whirl can be reduced with lubricating the rotor joints 
and splines. The whirl can be avoided with tight but 
nonlubricated interference fits. Joint locking with 
a high rotor unbalance is not recommended. 

Proportional Whirl Forces 

Increasing whirl forces can be most destructive. 
The simplest representation of the forces is a propor- 
tional relationship which couples rotor deflections by 
skewed-symmetric matrices to transverse whirl forces. 
Rotor hysteresis, 6 seal gaps,-^ journal bearings, ^ tur- 
bine blade bypass,® and impeller interactions drive 
usually subsynchronous forward whirls. External rotor 
rubbing drives a backward whirl that for hard rubbing 
can become supersynchronous (Fig. 1). 

Rotor Stabilization 


The purpose of rotor stabilization is to negate 
the whirl regeneration of proportional forces while 
reducing the whirl limit cycle caused by constant 
forces. Rotor suspension damping, stiffness, and 
asymmetry are whirl stabilizers ordered after their 
efficiency. Natural damping is usually too low to 
prevent whirl in high pressure applications. There- 
fore, squeeze film dampers are used in turbomachinery 
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not only for vibration attenuation but also for whirl 
prevention.^ However, such dampers are unsuitable for 
cryogenic applications. The difficulties prompted the 
development of damping seals. ® 

Damping seals are like annular seals, except for 
rough inner surfaces on stators (Fig. 2). Annular 
seals are a first step and were implemented to elimin- 
ate whirl from the high pressure turbopump (HPFTP) of 
the Space Shuttle Main Engine (SSME).2,3 However, 
such seals can destabilize a rotor resonance below 1/2 
the speed frequency. Damping seals improve this ratio 
to < 1/4, reduce leakage, and increase damping. Damp- 
ing seals were recently installed in a high pressure 
oxidizer turbopump (HPOTP) for SSME development as 
shown in Figure 3. 

Stability Criterion 

The stability condition of a single mass rotor 
(Fig. 4) suffices to examine the seal effect. Also, 
more complex models are represented, considering one 
resonant mode that is only negligibly detuned by the 
seal. Equation (1) describes the internal seal whirl 
effect with the whirl factor Q proportional to the 
seal damping C. The Couette factor c of equation (2) 
is the ratio of circumferential bulk flow velocity v 
to rotor surface velocity w. ft is the rotor speed in 
rad/s. Equation (3) is the characteristic equation of 
the model in Figure 4. Conformal mapping, 10 substitut- 
ing with equation (1), and reordering produce the 
inequality (4) which sets the speed limit. Ky and K z 
are the lateral suspension stiffnesses in y and z, m s 
is the rotor mass, and s is the Laplace operator, fti 
in equation (5) defines the first critical speed or 
resonance . 

Q = cftC (1) 


c = v/w (2) 

K K +Q 2 + sC (K +K ) + s 2 [C 2 + m (K +K ) ] 
yz x y z s y z 

+ s^ 2Cm + s^m 2 = 0 (3) 

s s 
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Figure 3. HPOTP Damping Seal Retrofit. 
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Figure 4. Simple Rotor Model. 
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( 4 ) 


« < [(K y +K z )/2m s + (K y -K z ) 2 /4C 2 ] 1 / 2 / c 
a i = [(K y +K z )/2m s ] 1/2 (5) 

The stability limit of inequality (4) is high for 
low c values, high rotor resonances, and asymmetric 
rotor suspensions. Asymmetry can be coincidental 
which makes symmetry a safer assumption for predicting 
speed limits. External whirl effects, such as turbine 
blade bypass, can be lumped into Q by raising c. Con- 
versely, external damping, such as from bearings, can 
be lumped into C by reducing c to maintain the same Q. 
Thus external effects are equivalently represented 
internally by the Couette factor c which becomes a 
measure of stability for comparing seals and other 
rotor dynamic components. 

Flow Friction 

The flow velocities relative to the rotor and the 
stator differ because of rotation. The velocities 
describe bulk flow so turbulent flow can be treated 
similar to pipe flow. Equations (6) and (7) describe 
shear stresses r r and x s respective to the rotor and 
stator. f r and f s are respectively the rotor and 
stator surface frictions which generally differ, u is 
the axial bulk flow velocity and p is the fluid mass 
density. The seal gap pressure gradients are given in 
equations (8) and (9) for the axial and circumferential 
directions x and q respectively, h is the gap thick- 
ness, p the pressure, and t the time. Note the .pro- 
jections of the shear stresses from equations (6) and 
(7) . The Reynolds numbers R r and R g in equations (10) 
and (11) are taken in the respective bulk flow direc- 
tions of the rotor and the stator, y is the fluid 
viscosity. The respective rotor and stator friction 
factors £ r and f s of equations (12) and (13) are 
Moody’s approximation with 2h as hydraulically equiva- 
lent pipe diameter D.H The effective roughnesses k r 
and k s reach pocket depth values for well circulating 
pockets. The triangular pockets of Figure 2 should be 
best for obtaining high k s values. 
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( 6 ) 


T r = p[u 2 +(v-w) 2 ]f r /8 


T = p (u 2 + v 2 )f /8 
s s 


(7) 


-3p/3x = p(uf (u 2 +(v-w) 2 )^ 2 + uf (u 2 + v 2 )^ 2 ]8h 
J- s 


+ p du/dt (8) 

-3p/3q = p[(v-w)f (u 2 +(v-w) 2 ) + vf (u 2 + v 2 )^ 2 ]/8h 

i. s 

+ p dv/dt (9) 

R r = p2h[u 2 +(v-w) 2 ] 1/2 /y (10) 

R g = p2h(u 2 + v 2 ) 1/,2 /y (11) 

f = 0.0055 + 0.55(k r /100h + 1/R r ) 1/3 (12) 

f = 0.0055 + 0.55(k g /100h + 1/R S ) 1/3 (13) 

Couette Flow 


The rotor generates a Couette flow which is 
established downstream of the seal inlet according to 
equation (14) which follows from equation (9) for the 
case of negligible circumferential pressure gradients 
and flow accelerations. Equations (15) and (16) are 
the asymptodes for u = 0 and u = 00 which enclose the 
Couette factor c for all surface friction ratios. 
Figure 5 shows c as a function of surface friction 
ratio and the axial velocity ratio u/w. c is low for 
high stator friction. The friction ratio is approxi- 
mately the cubic root of the roughness ratio, as the 
equations (12) and (13) indicate. Therefore, large 
roughness ratios are required. 

(v-w)f [u 2 +(v-w) 2 ] 1/2 + vf (u 2 +v 2 ) 1/2 = 0 (14) 

r s 

u = 0: c = l/[l+(f /f ) 1/2 ] (15) 

S L 
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Figure 5. 


Couette factor c = v/w versus the 


friction and axial velocity ratios. 



The integral dynamic effect of a seal is described 

fe 3 r stiffness K, damping C s mass m, whirl coupling 
cfiC* , and gyro coupling cfim’ as shown in vector equa- 
tion ( 17 ) » 



K+sC+s 2 m cft(C'+sm') 

-cflCC'+sm*) K+sC+s 2 m 



(17) 


The detailed derivations are given in reference 9, 

except for minor differences in rearranging frequency 
square terms herein. Some parts of the derivation are 
repeated for continuity. Equations (18) through (22) 
define the dynamic parameters and their response to 
the whirl frequency m. The time constants. are small. 
The basic . stiffness K 0 , damping C Q3 and induced mass 
xa 0 are given in equations (23) through (25)’ for inlet 
loss f]_; resistivities F s G s and H; dimensions h s L, 
and r of seal gap , seal lengths and seal radius; axial 
velocity u; and mass density p. The resistivities are 
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given by the equations (26) through (32) . The time 
constants are defined by the equations (33) through 
(36). 


K = [K o -(cf2) 2 m o ]/(l-ho 2 x 2 ) 

(18) 

C = [C (1+oj 2 tt )-kn to 2 (x-T, )-tK ]/(1+u 2 t 2 ) 
O 3.0 DO 

(19) 

m = [m (1+xT, (u^+^n) 2 )) - C (t-t ) ] / (l+w 2 x 2 ) 

0 D 0 3- 

(20) 

C' = [C (1+u 2 tt, )+2m w 2 (t-t, ) ] / (H-w 2 t 2 ) 
O DO D 

(21) 

m* = [2m (1 +oi 2 tt, )-C (t-t ) ] / (l-fu> 2 x 2 ) 

0 D O 3 

(22) 

K = (irr /h) u 2 ( 1+f . ) xL (F+H) / 2 
0 1 

(23) 

C = (irr/h) (t/t )L 3 (F4€)/12 

0 3 

(24) 

m = C t x /t, 
o o o a b 

(25) 

F = p [ f (u 2 +(v-w) 2 ) 1/2 + f (u 2 +v 2 ) i/2 ]/8h 

1 T S 

G = pu 2 [(f -g )/(u 2 +(v-w) 2 ) 1/2 
r r 

(26) 

+ (f -g J/(u 2 +v 2 ) i/2 j/8h 

S o 

(27) 

H = p[h (u 2 +(v-w) 2 ) 1/2 + h (u 2 -W 2 ) 1/2 ]/8h 

IT S 

(28) 

g r = h / (1+R k /100h) 

(29) 

g s « h s /(l+R s k s /100h) 

(30) 

h = (0.55/3) “(k /100h + 1/R ) 1/3 

X xx 

(31) 

h = (0.55/3)* (k /100h + 1/R ) 1/3 
s s s 

(32) 
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(33) 


x = pL/[pu(l+f 1 >+L(F+G) ] 


= pL/[pu(l+f 1 ) (4+6pu/L (F+G)+L (F+G) ] 

(34) 

= pL/ [4pu(l+f L )+L(F+G) ] 

(35) 

= p / (F+G) 

(36) 


The flow induced masses in equation (17) are 
dipol masses that are negligible against the rotor 
mass. The whirl coupling is indeed proportional to 
the Couette factor c as shown earlier in equation (1) . 
Table 1 gives the dimensions of diverse inlet and out- 
let seals of the HPOTP preburner pump impeller (Fig. 

3 ) . Table 2 lists the analytical results for laby- 
rinth, damping, and smooth annular seals. Table 2 
shows that the leakage rates of labyrinth and damping 
seals are comparable, while the smooth annular seals 
leak 1.6 times more. Tests by Rocketdyne produced 
the same ratio. The whirl ratio of damping seals is 
< 1/4 but 0.5 for the other seals. Damping seals have 
the highest damping while the smooth annular seals 
show the highest stiffness. 

Table 1. HPOTP Preburner Pump Seal Dimensions for 
Inlet and Outlet at Full Power 30,367 RPM 


PBP Seals 

Gap 

Diam. 

Leng. 

Pocket 

Press 

Dimension 

milsl 

inch 

inch 

mils 2 

psi 3 

Inlet Laby. 

6 

2.85 

0.63 

230 

2006 

Inlet Damp. 

6 

2.898 

0.66 

10(20) 

2006 

Inlet Smooth 

6 

2.898 

0.66 

0.032 

2006 

Outlet Laby. 

4 

3.345 

0.87 

80 

6505 

Outlet Damp. 

4 

3.3371 

0.88 

10(20) 

6505 

Outlet Smooth 

4 

3.3371 

0.88 

0.032 

6505 


NOTE: 1 = radial gap, 2 - effective pocket depth 

(actual depth), 3 = pressure difference. 
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Table 2. HPOTP Preburner Pump Seals 5 Dynamic for 
Inlet and Outlet at Full Power 30,367 RPM 


PBP Seals 

Stif . 

Whirl 

Damp . 

Ratio 

Leak 

Dimension 

Klb/in 

Klb/in 

Lbs /in 

c % 

Lb/ s 

Inlet Laby. 

2 

2 

1 

50 

4.94 

Inlet Damp. 

106 

46 

57 

23 

5.18 

Inlet Smooth 

151 

105 

53 

50 

7.90 

Outlet Laby. 

20 

36 

23 

50 

5.13 

Outlet Damp. 

447 

281 

388 

22 

5.09 

Outlet Smooth 

935 

544 

302 

50 

8.51 


NOTE: % ratio times maximum speed must be below 1st 

critical speed. Analysis is for 0.25 inlet 
loss, 2.59E-8 psi-s lox viscosity, and 
1.05E-4 psi(s/in.) 2 lox density. 

Seal Installation 

Figure 3 illustrates a retrofitting case where 
inlet and outlet labyrinth seals were replaced by 
damping seals to inhibit a 450 Hz whirl at 30,000 
RPM. Damping seals fit always into the space of 
labyrinth seals. Leakages are comparable, but stiff- 
ness and damping are much higher as seen in Table 2. 
Both changes inhibit whirl, lower vibrations, and 
increase critical speed. Damping seals seem to be 
more effective than squeeze film dampers, simpler, and 
smaller. Damping seals are suggested instead of 
smooth annular seals to reduce leaks and to suppress 
whirl. 

The high stiffness of damping seals must be, 
however, considered when replacing labyrinth seals. 
Eccentricities between bearings and seals can cause 
unexpected side loads. 

Soft Seals 


The side loads can be reduced with small eccen- 
tricities, if tight tolerance costs are not prohibi- 
tive. Soft seals are another concept worth 
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considering. Damping seals become softer as rough- 
ness is increased. Axial flow is thus reduced and the 
fluid compressibility much enhanced. Another way is 
the insertion of a labyrinth seal downstream to shorten 
a damping seal without changing leakage. 

Hard Seals 

The hard seal approach takes advantage of high 
seal stiffness and uses the seals for supporting the 
rotor. Critical speeds are so increased and loads of 
pump and turbine discharges are shared by the seals. 
Tight tolerances (eccentricities) can be circumvented 
by loose fittings for the bearings' outer races. The 
bearings support the rotor at startup and cutoff when 
the seals are unpressurized or soft. At operation 
speed the seals are stiff and can take over. Good 
bearing life is expected because the bearings float 
mostly during their operation in a loadless condition. 
Outer race spinning due to floating can be prevented 
with antirotation devices. 

Conclusion 

Damping seals are new technology for improving 
turbomachinery operations. The shown properties 
encourage many applications from retrofitting machines 
that are beset by whirl problems to new and more 
efficient designs. We are grateful and do not take 
this progress for granted - be blessed the same way. 
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